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Abstract. The complex Ginzburg-Landau equation is a universal model which has been 
used extensively to describe vaiious non-equilibrium phenomena. In the context of lasers, it 
models the dynamics of a pulse by averaging over the effects that take place inside the cavity. 
Ti:sapphire femtosecond lasers, however, produce pulses that undergo significant changes in 
different parts of the cavity during each round-trip. The dynamics of such pulses is therefore 
not adequately described by an average model that does not take such changes into account. 
The purpose of this work is severalfold. First we introduce the dispersion-managed Ginzburg- 
Landau equation (DMGLE) as an average model that describes the long-term dynamics of 
systems chai'acterized by rapid variations of dispersion, nonlinearity and gain in a general 
setting, and we study the properties of the equation. We then explain how in particular 
the DMGLE arises for Ti:sapphire femtosecond lasers and we characterize its solutions. In 
particular, we show that, for moderate values of the gain/loss parameters, the solutions of the 
DMGLE are well approximated by those of the dispersion-managed nonlinear Schrodinger 
equation (DMNLSE), and the main effect of gain and loss dynamics is simply to select one 
among the one-parameter family of solutions of the DMNLSE. 
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1. Introduction 

The complex Ginzburg-Landau equation (CGLE) is a universal model that governs the non- 
equilibrium dynamics of weakly nonlinear systems in the presence of gain, saturation, as well 
as linear and nonlinear dispersion, and describes a large variety of physical phenomena, from 
nonlinear waves to second-order phase transitions, superconductivity superfluidity, Bose- 
Einstein condensation, liquid crystals and strings in field theory [14, 21]. In particular, 
equations of CGLE type have often been used as a model for the description of mode-locked 
lasers [27, 28, 36, 39, 40, 52, 59]. The most well-known example in this context is referred to 
as the master equation of passive mode-locking [27, 36]. 

Lasers have a number of important scientific and technological applications, of course. 
In particular, the latest generation of Ti:sapphire femtosecond lasers [22, 65] could be used 
in spectroscopy, frequency metrology, and optical atomic clocks [23, 31, 47, 62]. One of 
the features of these lasers is the presence of large variations of dispersion and nonlinearity 
inside the cavity [19]. The two phenomena are called respectively dispersion management and 
nonlinearity management by analogy with optical fiber communications. Because of these 
variations, another mathematical tool which has been recently used in studies of femtosecond 
lasers is the dispersion-managed nonlinear Schrodinger equation (DMNLSE), which is a 
universal model that describes the long-term dynamics of weakly nonlinear dispersive pulses 
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subject to large, periodic variations of the nonlinear and/or dispersion coefficients [1, 24]. 
It should be clear, however, that the CGLE and DMNLSE are both inadequate descriptions 
of Ti:sapphire lasers. On one hand, the CGLE is derived under the assumption that the pulse 
changes per round-trip inside the cavity are small, a condition which is grossly violated in 
femtosecond lasers. As a result, the CGLE cannot provide an accurate description of the 
dynamics of such lasers. On the other hand, gain/loss effects are obviously significant in any 
laser, so the DMNLSE canmiot be a fully accurate description of a laser. 

In this work we address both of the above issues. Starting from a multi-dimensional 
perturbed CGLE with time-dependent coefficients, we derive a new equation — which we 
call the dispersion-managed Ginzburg-Landau equation (DMGLE) — that describes the long- 
term dynamics of weakly nonlinear systems subject to rapidly varying dispersion, nonUnearity 
and gain/loss. We then study the properties of this equation and we use it to investigate the 
behavior of dispersion-managed soUtons in femtosecond lasers. The outline of this work 
is the following: In section 2 we obtain the DMGLE in a general setting as the equation 
which describes the average dynamics of systems governed by the CGLE with large and 
rapidly varying coefficients. In section 3 we then show how in particular the DMGLE arises 
in the description of Ti: sapphire femtosecond lasers. In section 4 we study the symmetries, 
rate equations and linear modes of the DMGLE, and in section 5 we discuss solitary-wave 
solutions and their properties. Section 6 concludes with some final remarks. 

2. From the CGLE to the DMGLE 

Here we obtain the DMGLE in a general setting as the equation which describes the long-time 
dynamics of multidimensional systems governed by the CGLE with large, time-dependent and 
rapidly varying coefficients. We start from a perturbed cubic CGLE in dimensionless form 
with time-dependent coefficients 

\p^^/t'',t)Vh + n{t/ta,t)\q\^q = ig{t/ta,t)q, (2.1) 

where q = q{x,t) € C, x = (xi , . . . ,xn) g M'^, f e K, is the Laplacian operator, and the 
coefficients p and n are complex: p{-) = pre( ) + 'Pim( ) etc. In the context of mode-locked 
lasers, these quantities have the following physical interpretation: 

• Prc{t /ta,t) quantifies chromatic dispersion, 

• Pim {t/ta,t) quantifies spectral filtering or band-hmited gain, 

• nre{t/ta,t) quantifies nonUnear interactions mediated by the Kerr effect, 

• n\m{t/ta,t) quantifies gain saturation, 

• gvc(t/ta,t) quantifies linear gain 

(see section 3 for a more detailed discussion). Here is a short temporal scale characteristic 
of the problem, which will be defined more precisely below. Without loss of generality, the 
function g(t/ta) can be taken to be real, since an imaginary component of g{t/ta) can always 
be eliminated by defining a rescaled field via the phase exp [;/ ^im(?/?a)d/] . But the same 
statement does not apply to p{t/ta) and n{t /to), of course. 

Solutions of (2.1) with complex coefficients have been extensively studied; e.g., see 
Refs. [12, 14] and references therein. These studies treat the standard CGLE, namely, the 
case in which all the coefficients are constant. Solutions of (2.1) with large, time-dependent 
and rapidly varing coefficients, on the other hand, can be very different from those of the 
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standard CGLE, and characterizing their behavior is a nontrivial task. Moreover, to the best of 
our knowledge, such problem has not been addressed in the literature, despite its theoretical 
and practical relevance. Our goal is precisely to elucidate the dynamics of such solutions. 
We do so by employing a multiple- scale analysis similar to that used to study optical fiber 
communication systems in the presence of dispersion management [1, 63]. 

Multiple scales expansion and the DMGLE. If the temporal variations of the coefficients 
in (2.1) are not too large, one can employ a standard multiple-scale perturbation expansion, as 
was done for the NLSE in [63]. The result is that, to leading-order, the pulse satisfies a CGLE 
with constant coefficients that are simply the averages of those appearing in (2.1), Higher- 
order terms then give corrections to this leading-order behavior. (Of course different starting 
models can give rise to different scenarios, e.g., see Refs. [45, 55].) When the coefficients 
in (2.1) exhibit large temporal variations, however, the situation is not as simple, since the 
standard perturbation expansion breaks down, as it does for the NLSE [63]. It is this latter 
case that is of interest here. There are three key assumptions for our analysis: 

(i) All of the coefficients in (2.1) are periodic with respect to their first argument. That is, 
p{ll,+\,t) = piC,t) VC e R and similarly forn( • ) andg(-). 

(ii) The coefficients in (2.1) are rapidly varying. That is, the dimensionless period is ?a <C 1. 

(iii) The function p{-) can be decomposed into the sum of two terms, describing respectively 
large periodic zero-mean oscillations and an 0(1) residual. Namely, 

p{t/ta,t) = ^Ap{t/ta)+Pies{t/ta,t), (2.2a) 
ta 

with 

'jAp{t/ta)dt = 0. {2.2b) 



Note that both and pres can have both a real and an imaginary component. (Indeed, this 
will be the case in section 3.) Note also that the double time dependence of p{-), «(•) and g{-) 
in (2.1) allows one to capture cases in which the corresponding physical parameters have both 
a fast and a slowly varying component with respect to time. Even though the simplest case is 
when all of the coefficients average to a constant (as in optical fiber communications [1]), the 
more general situation does not introduce significant complications in the method. 

With the above assumptions, we introduce the fast temporal scale = t/ta and look for 
solutions in the form 

^(x, C,f) = C,f) +f«<7(')(x, CO + . (2.3) 

It is convenient to denote by f{x,t) the average of a function f{x,t/ta,t) over one map period, 
namely: 

/(x,0 = //(x,C,OdC = yjf{x,t'/ta,t)dt'. (2.4) 

'a 

We then substitute the ansatz (2.3) into the CGLE (2.1) and solve the resulting perturbation 
expansion following standard procedures. At leading order in the expansion we have the 
following linear problem: 



(2.5) 
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The solution of (2.5) is of course trivially obtained using Fourier transforms. We denote 
Fourier transforms by a circumflex accent throughout, with the transform pair defined as 

/(k) = j[/(x)] = /e*-/(x) (dx) , (2.6a) 

/(X) = J- V»] = ^ /e-*-''/(k) (dk) , (2.6b) 

where (dx) = dxj • --Axn is the volume element in M!^ . (Integrals are complete — meaning 
over all of — throughout this work unless otherwise noted.) With these conventions, the 
general solution of (2.5) is 

^W(k,C,0=Z(k,S(O)M(k,0, z(k,^)=e-'«'^'/2, (2.7) 

where A:^ = k • k, 

S(C)=S, + /A/,(ndC' (2.8) 



and S„ is an arbitrary integration constant. Equation (2.7) separates the fast dynamics of 
the pulse from the slow dynamics: The exponential factor in (2.7) accounts for the rapid 
"breathing" (periodic compression/expansion) of the pulse, while the slowly varying envelope 
M(k,f) encodes the information about the core pulse shape. Inverting the Fourier transforms 
we get the leading-order solution in the spatial domain, 

1x2/2^ 

5(«)(x, CO = l(x,S(0) * «(x,0 , ;t(x,^ ) = , (2.9) 

V(27r0"<^ 

where = x • x and y/z is taken on the principal branch (with | arg(z)| < n), and where the 
asterisk denotes the convolution integral, defined as 

= Ifix )g(x - x') (dx') , (2.10) 

Above and below we used the convolution theorem, stating that ^ ' [fg] =f*g and ^ [fg] = 
f*g /{InY ■ Note that ^(k,0) = 1 and ;if (x,0) = 5(x), where 5(x) is the A^-dimensional 
Dirac delta. Therefore, at those points for which E(Q = 0, the solution q^^\x,l^,t) 
coincides with u{x,t). Note also that, since is even in x, the parity of ^(^^(x,-,-) 

is determined by that of m(x, •). That is, if u{x,t) is even with respect to x, ^("^ (x, ^,t) is also 
even; vice versa, if u(x,t) is odd ^^"^(x, ^,t) is also odd. 

The function u{x,t) is arbitrary at this stage, and as usual must be determined at higher 
order in the expansion. At (9(1) we have the following forced linear problem: 

'^■^''+^Aes(C,0v,2^(«)+n(C,0k'''¥^''')-<-^(C,0^(°' 

(2.11) 

The solution of (2.11) can also be easily obtained using Fourier transforms. To avoid 
secularities, however, a Fredholm solvability condition must be satisfied (which as usual 
expresses the orthogonahty of the forcing to the solutions of the homogeneous problem). 
After tedious but straightforward algebra, this condition yields 

+ ^P(0'^x^« + //«(x+x')«(x+x")«(x+x'+x")^(x'-x") (dx')(dx") = igit)u, 

(2.12) 
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or, equivalently, in the Fourier domain, 

~ ^^(0^« + //«(k+k')«(k+k")"(k+k'+k")%'-k") (dk')(dk") = ig{t)u., 

(2.13) 

where for brevity we have introduced the shorthand notation mj^) = u{x,t) and Mf^) — M(k,f). 
(That is, subscripts in parenthesis denote functional dependence, not partial differentiation.) 
The superscript * denotes complex conjugation throughout, and the integration kernels (which 
are in general complex and time-dependent) are given by 

m,t) = ^^/n(C,Oe'^(^)«dC (2.14a) 

and 

K{x'-x",t) = (27r)2V"'[^(k'-k",f)]- (2.14/7) 

[Again, a straightforward calculation shows that, as in the one-dimensional DMNLSE [1], 
K{-,t) depends only on the dot product x' • x", not on x' and x" separately. Note also 
that p[t) = Ptc&{t)-\ We refer to (2.12) [or equivalently (2.13)] as the dispersion-managed 
Ginzburg-Landau equation, or DMGLE. 

Remarks. Several comments are now in order: 

• If A/?(Q = 0, then choosing E„ = one has E(Q =0, which implies ^(i^,f) = l/(27r)^'^ 
and k{x! ■ x",t) = 5{x' - x") . Then (2.12) and (2.13) reduce to the usual CGLE with 
constant coefficients written in the spatial and Fourier domains. 

• If Pim{^,t) =nim{^,t) =0, equation (2. 12) reduces to the multidimensional extension of 
DMNLSE [1, 24], perturbed by the addition of linear gain and loss when g{t) ^ 0. Then, 
if Ap{Q = 0, equation (2.12) further reduces to the NLSE with constant coefficients, 
again possibly with gain and loss. 

• Since the NLSE, the CGLE and the DMNLSE arise in many different physical situations, 
one can expect that the DMGLE will also have a wide range of applicabihty. Indeed, 
the above formulation is fairly general, and it allows one to study a variety of different 
scenarios, including spectral filtering and nonlinear loss and gain. 

• Equations of NLS type or CGLE type with nonlocal nonlinearities arise in various 
physical situations. In particular, apart from optical fiber communications, they also 
appear in water waves, in which case they are known to be integrable [53]. 

• Like the DMNLSE (its counterpart for conservative systems), the DMGLE is a reduced 
model that retains the essential features of dispersion-managed systems while bypassing 
the complicated dynamics that take place within each period. As such, it should prove to 
be a useful model to investigate the long-time behavior of nonconservative dispersion- 
managed systems, like the DMNLSE did for conservative ones. 

• Since q{x, ^^,,t) = u{x,t) whenever = 0, the value of Zo simply dictates at which 
points inside the map the full solution of (2.1) coincides with its core pulse shape. At the 
same time, however, the choice of S„ appears in , t) and its transform, and therefore 
can sometimes be exploited to make them simpler. (For example, the corresponding 
kernels for the DMNLSE in the lossless case can be made real [1].) 

• If 7^ const, the kernels K{^,t) and K{^,t) can have a nonzero imaginary part even 
when Apini(^,f) = 0. Indeed, this is the case in optical fiber communication systems 
with loss compensated by periodically spaced amplifiers, when one takes into account 
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that the effective nonlinear coefficient in the NLSE is distance-dependent due to the 
periodic power variations originated by the loss/amplification cycle. 

• A multiple-scale averaging similar to the one described above can also be performed 
for the quintic Ginzburg-Landau equation. The result of such an averaging is that the 
quintic term in the original equation produces a quintic term in the DMGLE with a 
four-dimensional convolution integral and a corresponding integration kernel. Such a 
situation will not be considered here, however. 

Finally, we reiterate that A/7(Q need not be real for the purposes of the multiple- scale 
expansion: the derivation applies without change when Apini(C) 7^ 0- Of course, some 
of the conclusions are different: S(C) acquires an imaginary part, and this has two main 
consequences: (i) it leads to a periodic ampUfication/suppression of high-wavenumber 
components in the fast dynamics; (ii) it leads to different kernels in the DMGLE. While 
the first of these consequences can be considered to be minor (since it has no effect on the 
long-term behavior of solutions), the second can in principle produce different and interesting 
phenomena. 

3. The DMGLE in femtosecond lasers 

We now discuss how the DMGLE appears in the description of Ti:sapphire lasers. Figure 1 
shows a diagram of a typical experimental setup. The system consists of a continuous-wave 
(CW) pump, a Titanium-doped sapphire crystal and a set of prisms and/or mirrors. The 
crystal, which constitutes the nonlinear medium, is characterized by a Kerr-type nonhnear 
response as well as by gain and gain saturation, and has a large normal group velocity 
dispersion. [That is, for CW waves whose electric field is written as Re(£' e'*^"'^™), where k{(o) 
represents the linear dispersion relation it is k!' = d^k/doP' < 0.] The prisms and mirrors, on 
the other hand, are especially designed to have a large anomalous dispersion to compensate 
that of the crystal. The output of such a system is a periodic stream of optical pulses with a 
typical duration of about 10 fs, a spectrum typically centered at 830 nm and 70 nm wide, and 
with a typical repetition rate of 90 MHz (e.g., see Ref. [31]). 

Laser model. The evolution of a quasi-monochromatic optical pulse in a Kerr-type medium 
is usually well described by the NLSE, possibly with varying coefficients [10, 19, 57]. To 
accurately capture the behavior of pulses in a laser, however, it is obviously necessary take 
into account the gain and loss dynamics [19, 46|. In Ti: sapphire lasers, the gain dynamics 
(including hnear gain and gain saturation) takes place inside the crystal, while the losses are 
concentrated in the reflecting mirror(s) and the output coupler. Since the precise dependence 
of the gain on the amplitude is compUcated and in general not known in closed form, one is 
forced to approximate. This is typically done by characterizing the response of the nonlinear 
medium via two main parameters: the smaU signal gain, Gss (namely, the gain experienced 
by signals of small ampUtude) and the saturated gain, G"sat (namely, the gain experienced by 
signals of large amplitude). Of these two parameters, Gsat can be obtained indirectly, because, 
for each system configuration, it must compensate exactly the total losses per cavity round- 
trip. The value of Gss, however, is hard to characterize theoretically or experimentally. Note 
that the value of Gsat depends on the specific system configuration, most notably on the pump 
power. Note also that gain relaxation is also present, but occurs over much slower temporal 
scales than those of a single optical pulse. Since we are interested in the behavior of pulses 
near steady-state, this difference in the time scales is not a concern. 
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Figure 1. Schematic diagram of a prototypical Ti:sapphire femtosecond laser. See text for a 
detailed description. 



The gain dynamics is typically approximated by choosing a heuristic gain response 
function that interpolates between the small signal gain at low powers and the saturated gain 
at high powers. In particular, it was shown that choosing a Lorentzian-shape function yields 
a good approximation for the gain saturation (e.g., see Ref. [19, 27]). Let E{Z,T) be the 
slowly varying envelope of the electric field, Z = ziab is the physical propagation distance and 
T = ?iab — ziab / Vg is the retarted time [that is, the time in a reference frame that moves with 
the group velocity Vg = 1/A:'(a)) of the pulse]. Neglecting all other effects for the moment, the 
combination of linear gain and gain saturation yields the following equation for the evolution 
of the electric field envelope inside the Ti:sapphire crystal: 



'dZ 



1 + \E/Eo\^ 



(3.1) 



where Eo represents some appropriate reference amplitude. As in [35], we now expand the 
fraction 1/(1 + \E/Eo\^) in Taylor series near E — 0. Retaining terms up to next-to-leading 
order, the right-hand side of (3.1) is then replaced by Gss(l — |£^P/|£^oP)£^- 

Assuming that the gain spectrum does not change appreciably with the pump power over 
the operating range of the laser (which is a good approximation in Ti: sapphire lasers), one can 
neglect the second derivative of the cubic term when reinstating the time dependence. Then, 
recalling that the gain dynamics is confined to the crystal and combining the gain dynamics 
with the dispersive and nonlinear effects and the loss arising of the mirrors yields the following 
CGLE in dimensional form with distance-dependent coefficients: 



dE 1 



d^E 



^{k"{z)-im)-^^2 



(7(Z) +/Gni(Z))|£|2£ = /G„(Z)£, (3.2) 



where k"{Z) and 7(Z) are respectively the dispersion coefficient and the nonlinear coefficient, 
while j3 (Z) quantifies band-limited gain. 

The coefficients in (3.2) are conveneniently parametrized by introducing the "indicator" 
function /(Z), which equals 1 for values of Z corresponding to locations inside the Ti:sapphire 
crystal and otherwise. The coefficients in (3.2) can then be written as: 

k"{Z)^l^^,J{Z)^k'^,^{\~I{Z)), (3.3a) 

P (Z) = /3Ti:s/(Z) , 7(Z) = yri:s/(Z) , 0.3b) 



Go (Z) = Gss /(Z) + Gdcm ( 1 - /(Z) ) , Gnl (Z) = 



\Eo\ 



(3.3c) 



where the constants k'L.. 



^Ti s'^i:s > are the dispersion and nonlinear coefficient of the 
Ti:sapphire crystal; k"^^^ < quantifies the dispersion of the mirrors, prisms and output 
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coupler, as well as that of the air (which is experienced during the free-space propagation 
of the pulse in the cavity), and Gdcm < does the same for the loss. Note that, as written, 
(3.3a) and (3.3c) model the dispersion compensation and linear loss as taking place uniformly 
over the whole length of the cavity not occupied by the Ti:sapphire crystal. This is obviously 
not the case. Nonetheless, for the purposes of the averaging it is immaterial whether these 
effects occur at discrete locations or whether they are distributed, because all the processes 
occurring outside the crystal are linear. 

Nondimensionalization and the DMGLE. Next we proceed to nondimensionaUze (3.2). This 
requires more care than is necessary for the NLSE or the CGLE, because the dispersive effects 
are distance-dependent. It also requires more care than is necessary for the DMNLSE, because 
the nonlinear effects are also distance-dependent here. 

Let X = (flab — ziab / Vg) /f* and f — ziab / z* be respectively the dimensionless retarded time 
and the dimensionless propagation distance, where z« is a typical distance (defined below) and 
f* is a typical time scale (e.g., the pulse duration). Also, \e\.q = E/E^, be the dimensionless 
slowly-varying complex envelope of the electric field, where = \E^, p is a typical pulse 
power. Taking into account the fact that nonlinearity only acts inside the crystal, we then 
set z* = (i/ri:s/^a)zNL. where Lxiis is the length of the Ti:sapphire crystal. La is the total 
length of the cavity and znl = 1/ {Yvi-.sP*) is the distance at which nonlinear effects become 
relevant. These changes of dependent and independent variables, transform (3.2) into the 
one-dimensional version of (2.1), with 

p{Q = d{Q-ib{Q, d{Q^-k"{Q/k':, b{Q = bI{Q, (3.4a) 

n{Q = {l+ic)IiQ/I, giQ=GoiQz. (3Ab) 

(with a slight abuse of notation), where = //fa as before, and where fc" = tl/z*,b~ fhiJ/ K 
andc= (Gss — Gsat)/(yri:s|£'oP)- [Note that /^ =La/z* and / = Lxiis/i^a-] Without repeating 
the perturbation expansion of section 2, we can then apply its results to conclude that, to 
leading order, the dimensionless slowly varying envelope of the electric field is given by 

q{x,t) = J-^ [(2(A:,/)e'^('/''')*'/2] (3 5) 

with 

s(0 = / [d{i;') - d- i{b{C) - b)] dC , (3.6) 



and where the slowly-varying core u{x,t) = ;F^^[u{k,t)] solves the following DMGLE: 

du I , J r-, d^u 
dt 2^ ' dx^ 

+ (1 + ic) JJ M(;t+y)M(;,+y/)M(^_^y_^y,)/?(yy/) 6x'6x" = igu , (3.7) 

with b,c,d,g e M, and 

R{x'x") = JJ e-*'^-*"^V(A:'A:") 6k' 6k" , (3.8a) 

with 

If p{Q e K [i.e., b{Q = 0], it is r*(§) = r(-^) and R*{^) = R{-^). When b = c = g = 0, 
(3.7) reduces to the DMNLSE derived in [1, 6]. 

The multiple-scale expansion is of course formally justified only when the expansion 
parameter ta is small. Nonetheless, the DMNLSE — obtained under the same assumptions by 
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neglecting the gain and loss dynamics — has been shown to provide a good qualitative (and 
in some cases even quantitative) description of the actual behavior of pulses in Ti:sapphire 
femtosecond lasers even for system configurations where ta = 0{l) [57]. One therefore 
expects that the DMGLE will also be a good model even in those situations. 

As in section 2, the coefficients b, c, J and | could also depend on the evolution variable t. 
Such situations occur when the average dispersion, gain and gain saturation exhibit variations 
over slower temporal scales compared to the characteristic duration of the pulses in the cavity. 
Note also that the average dispersion depends on the details of the Kerr-lens process in the 
crystal, on the mirrors, the prisms and on the total cavity length. If d = const, one can make 
its value unity by choosing t^ = k"z*. However, the general form is convenient if one wants to 
compare systems having different values of average dispersion (as in Ref . [57]), since in this 
case one can do so within the framework of a single DGMLE, without having to go back to 
the COLE and choose different normalizations for each case. 

Two-step maps. For the piecewise-constant two-step maps defined by (3 .3), the kernels r{^) 
and R{^ ) in (3.7) can be computed exphcitly, and assume a very simple form. Taking the 
origin of the map at the output coupler, we have 

/m-l' Ce [0, 0/2) u [1-0/2,1), 

\i [0/2,1-0/2), ^ ^ ^ 

where < < 1 represents the fraction of the cavity length not occupied by the Ti: sapphire 
crystal. This yields I = 1 — . Recalling that the linear effects acting on the pulse can 
be distributed throughout the portion of the cavity not occupied by the crystal, we then 
parametrize the zero- mean part of the dispersion as [1, 57] 

fV20 Ce [0,0/2) U [1-0/2,1), 

\-^v/[2(l-0)] Ce [0/2,1-0/2), 

where the map strength parameters quantifies the Li-norm of Ad(Q, as discussed in section 4. 
(This definition differs from the previous one [1, 57] by a factor 4.) Then 



's+C [0,0/2), 

S_(C-l/2) Ce [0/2,1-0/2), (3.11) 

z+iC-i) Ce [1-0/2,1), 



where S+ = 2{s + ibo)/9 and S_ = —2{s + ibo)/{l — 0), and with bo = bdta/2. After some 
straightforward calculations, one then obtains the kernels as 

1 

(2"^)2 



K^) = 7^ smc[{s + 2ibo)^/4] , (3.12a) 



where sinc(>') = (sin >')/>■ and ci(>') = /~cosfdf/f is the cosine integral. The functional 
form of both kernels is identical to that of the DMNLSE. (The kernel in Ref. [57] differs 
by (3.12a) by an overall multipUcative constant due to the different choice of normalizations.) 
In the DMGLE, however, the kernels r{^) and R{^) appear multiplied by the complex factor 
l+ic, which has important consequences on the equation (as we show in section 4) and on 
its solutions (as we show in section 5). Note that when b = both kernels are real, as for 



The DMGLE and its application to femtosecond lasers 



10 



the DMNLSE without linear gain and loss [1], and the specific details of the functions M.{t) 
and I{t) do not affect the DMGLE nor its kernels. For this class of systems the DMGLE is 
a universal equation, in that it provides a unified description for the dynamics of different 
systems (like the NLSE, CGLE and DMNLSE). 

4. Symmetries, rate equations and linear modes 

We now look at the properties of the DMGLE. For simplicity, in the remainder of this work 
we will restrict our attention to the form of the DMGLE that appears in femtosecond lasers, 
and to the case in which the coefficients b, c, d and g appearing in (3.7) are independent 
of time. We emphasize however that many of the results that follow can be extended in a 
straightforward way to the more general form of the DMGLE (2.12). 

Symmetries. The DMGLE (3.7) enjoys a number of symmetries. Let u{x,t) be any solution 
of (3.7). Then Ue{x,t) is also a solution of (3.7) Ve e M, where: 

(i) Phase invariance: U£{x,t) = e'^u{x,t). (The same symmetry applies for the NLSE, 
CGLE and DMNLSE.) 

(ii) Space translations: Ue{x,t) = u{x— e,t). (The same symmetry applies for the NLSE, 
CGLE and DMNLSE.) 

(iii) Time translations. Ui;(xj)=u{x,t — e). (The same symmetry applies for the DMNLSE. 
It also applies for the NLSE and CGLE in the case of constant coefficients. In 
general (2.1) is not invariant under time translations, however.) 

(iv) "Chirp": Ue{k,t) = e^'^IH{k,t) is a solution of the equation with r{B, ) replaced by 
re(i^) = e'^^r(i^). (The name comes from the fact that quadratic variations of the pulse 
phase are referred to as "chirp" in optics [10, 48].) This symmetry, which is the same 
as that for the DMNLSE [1], has no analogue for either the NLSE or the CGLE. 

(v) Galilean boosts [generalized Galilean boosts if <i(-) is time-dependent]: if = 0, 

ue{x,t) = e'["^-5^''"(')lM(x- em(0,0 (4.1) 
withOT(f) =/oJ"(f')df'. (The NLSE and CGLE with varying coefficients admit a similar 
invariance [34]. For the DMNLSE, as well as for the NLSE and CGLE with constant 
coefficients, it is simply m{t) = dt. For the DMGLE with b^O, (4.1) appUes only for 
those solutions that can be analytically continued off the real x-axis.) 

(vi) Generalized scaling symmetry. If u{x,t;g,s) is a solution of (2.12), so is 

Ue(x,t;g,s) =au{ax,a^t;g/a^,a^s) (4.2) 

with a = l+e and where s is the map strength parameter, as discussed below. (The 
parametrization of (4.2) is chosen so that the value e = reproduces the original 
solution, as in all other cases.) 

These invariances will allow us in section 5 to construct a two- or three-parameter family 
of solutions of the DMGLE (depending on whether b ^ or b = 0, respectively) from a 
single stationary solution, in a similar way as for the NLSE, DMNLSE and CGLE. As we 
show later, however, unhke for the NLSE and DMNLSE (and hke for the CGLE) the scaling 
invariance does not generate a one-parameter family of solutions. Note also that, as a result 
of the chirp symmetry, the constant Eg in (2.8) can be chosen arbitrarily, and does not affect 
the solution q{x,t) of the original problem (2.1), since it simply amounts to a redefinition of 
u{k,t). (Proper choice of Zg can be useful to make the kernels simpler, however.) 
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Map strength. In the statement of the scaling invariance above we used the obvious fact that 
the solutions of the DMGLE (2. 12) depend parametrically on g. One should realize, however, 
that, Uke for the DMNLSE, solutions of the DMGLE also depend on a parameter s called 
the map strength, which quantifies the size of the zero-mean dispersion variations. The map 
strength (which for two-step maps was introduced in section 3) can be formally defined for 
any map as the Li-norm of Ad (in the DMGLE as for the DMNLSE [41]), namely: 

i = M = /IMO|d?. (4.3) 



One can then obtain expUcitly the dependence of the kernels r{B,,t) and R{x',x",t) on s by 
writing Ad{Q [and consequently S(Q via (2.8)] in terms of normaUzed functions. That 
is, given any choice of map Ado{Q, we can define the normalized function Adref{Q = 
Ado{Q/\Ado\. One can then introduce a one-parameter family of dispersion functions 

Ads{C)=sAdr,f{C). (4.4) 
In this way one can study the behavior of solutions for different values of map strengths 
within the framework of the DMGLE, without needing to go back to (2.1). When necessary, 
we express this dependence explicitly by writing the solution of (2.12) as u{x,t;g,s), like we 
did in (4.2). Of course, in the limit s — > the DMGLE (3.7) reduces to the CGLE with constant 
coefficients. Using the map strength parameter, it is now easy to show that, if d,g and K are 
independent of /, the generahzed scahng invariance (4.2) holds. Of course this invariance 
reduces to those of the DMNLSE, CGLE and NLSE in the cases s = and/or b = c = g = 0. 

Rate equations. Like the CGLE, and unlike the NLSE and DMNLSE, the DMGLE is not a 
Hamiltonian system. Hence it is not possible to use Noether's theorem to derive conservation 
laws from the symmetries of the equation as done for the NLSE [33] . Nonetheless, by analogy 
with the NLSE/DMNLSE, it is still possible to associate each symmetry with a rate equation. 
For brevity we only list the first three such equations, and only in the special case b = 0: 

^f\u\'^6x = 2gf\u\'^6x-2cfffu*^-^u^x+x')U{x+x'')U*x+x'+x'')^{x^^^^ (4-5a) 

— Im / u*dxudx = 2^1m / u*dxudx — 2clm J (^dxU*) U(^x+x")'^{x+x")ii{x+x'+x")^{x'x")'^''^"'^ ^ 

(4.5fc) 

^^['^1 \dxU\^dx- ffful^)U(^x+x')U{x+x")ulx+x'+x")R{x'x'')<ix^ = 

g [dj \dxu\^dx-2 Iu*(^f(x+x')U(x+x") ulx+x'+x")^ix'x")^-'('(ix"(i-'(] 

+ c[dReJJJ {d^u*)u(^x+x')U{x+x'')U*(x+x'+x")^ix'x")(ix'dx''dx 

+ 2/ 1 JJu(^x+x')U{x+x'')U*x+x'+x'')f^ix'x'')'^<ix''\^(ix] , (4.5c) 

where for brevity we used the notation dxU — du(^x)/^x etc. Each of the above rate equations 
reduces to those of the NLSE [33, 60], CGLE and DMNLSE in the appropriate Hmits, of 
course. (It appears however that, even in the simpler version obtained for s — 0, that is for 
the CGLE, the third rate equation is not widely known; cf. Ref. [13].) In particular, for the 
NLSE and DMNLSE, the right-hand- side is zero, and the integrals on the left-hand side of 
equations (4.5), which are then the conserved quantities of the corresponding equation, are 
respectively the pulse energy, the momentum and the Hamiltonian. Because of the physical 
meaning associated with these integrals, the rate equations for the NLSE and DMNLSE 
have been extensively used to study the evolution of various pulse characteristics. The rate 
equations of the DMGLE should therefore prove to be similarly useful. 
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Linear modes. If u{x,t) is any solution of (3.7) and Ue{x,t) = u{x,t) +ew{x,t) is also a 
solution, then w{x,t) belongs to the nullspace of the linearized DMGLE operator about u{x,t). 
That is, solves the linearized DMGLE L[w,m] = 0, where (for ^ = 0) 

xr 1 1 jd^W ._ .-^ rr * r. , / , // 

L[W,U\ =l-^ + 2'^-5^ -^^^ + 2(1 +«c)//W(;,+y)M(^+y/)M(^_^y_^y,)/?(yy/)dydx' 

+ (1 + ic) 6x'6x" . 

(4.6) 

It then follows that, for each continuous invariance of the DMGLE, there exists a solution of 
the hnearized DMGLE in the form 



dus{x,t) 

w{x,t) = 



de 



(4.7) 

6=0 



similarly to the CGLE, NLSE and DMNLSE [41]. Then, by simply applying the invariances 
listed above, we obtain the following set of Unear modes and generalized linear modes: 

du du 
wi=m, W2 = — ^ , W3 = — ^ , W4 = ixu, (4.8a) 
ox ot 

du du _du 
W5 =u+x-^+2s-^ -2g-^ , (4.8fe) 

corresponding respectively to phase rotations, space translations, Galilean boosts and scaling 
transformations as given above. (Of course W4 only applies when h = 0. Note also that the 
chirp symmetry requires changing the kernels, so the above analysis does not apply.) The 
above linear modes and generalized modes satisfy the relations 

L[wi,m] =L[w2,u] =L[w3,u] = 0, (4.9a) 
L[w4.u] = iw2 , L[w5 . u] = 2w3 . (4.9^?) 

In the special case of traveling wave solutions, time translations are simply a composition 
of space translations and phase rotations. As a result, for traveling wave solutions ws can 
be expressed as a linear combination of w\ and W2, and is not an independent mode. This is 
consistent with the familiar result [60] that for the NLSE in A' spatial dimensions the linearized 
operator around a solitary wave solution has a zero eigenvalue of multiplicity 2N+2. (The 
result also applies to the DMNLSE [54].) 

The linearized DMGLE is of course useful to study the stability of solutions. As 
in the NLSE and DMNLSE, the linear modes of the DMGLE generate changes in the 
solution parameters, and can be used to quantify the change of these parameters under 
perturbations [29, 41]. Because of the presence of a nontrivial nullspace, secularities arise 
even if all eigenvalues of L have zero or negative real part, and as a result the change 
in the solution of the DMGLE does not in general remain bounded in time. As usual, 
secular terms are removed by taking the solution parameters to be slowly dependent on time 
and determining their evolution by projecting the perturbation onto the linear modes. The 
existence of a generaUzed mode corresponding to the generaUzed scaling invariance is then 
important because the coupling between amplitude and phase in the NLSE/DMNLSE is the 
mechanism whereby the variance of noise-induced phase perturbations grows cubically [26] 
(similarly to the well-known coupUng between frequency and timing jitter [25]). Even though 
the relation between amplitude and phase is broken in the DMGLE due to the presence of 
gain and gain saturation, the existence of a generalized mode associated to the phase suggests 
that the noise-induced phase variance could still grow cubically. 
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5. Solutions, stability and parameter dependence 

We now discuss special solutions of the DMGLE. Again, for concreteness, we restrict our 
consideration to the specific form of the DMGLE arising in femtosecond lasers, namely (3.7) 
with constant coefficients. 

Soliton solutions of the DMGLE. We start by looking for stationary solutions, that is, 
solutions in the form 

M(x,0=/We'^''/^ (5.1) 

[Recall that a family of traveling wave solutions can be obtained from (5.1) by applying the 
invariances of the DMGLE.] Substituting the Fourier transform of the above ansatz in (3.7) 
yields the following nonlinear integral equation 

[^k^ + ^{d-ib)k^ + ig]f = 

(1 + ic) !! f(k+k').f(k+k").f(k+k'+k")''(k'k") dk'dk" . (5.2) 

Although closed-form solutions are not available unless ''( •)= const, (5.2) can be efficiently 
integrated numerically (see Appendix). Note also that, similarly to the DMNLSE [41, 44], 
fast numerical methods can be used for the calculation of the double integral in Eqs. (3.7) 
and (5.2) and thus to solve numerically both the integral equation (5.2) and the DMGLE (3.7) 
itself (again, see Appendix). Hence, the computational complexity of the DMGLE is no less 
and no greater than that of the original, un-averaged equation (2.1). 

Of course, solutions of (5.2) are not sohtons in the mathematical sense of the term 
(that is, solutions corresponding to the discrete spectrum of the scattering problem associated 
to the given nonlinear partial differential equation [9]), but rather solitary waves. As common 
in physics and optics, however [10, 18, 36, 48, 65], we will refer to such a pulse as a soliton 
— or, in the case of the DMNLSE and DMGLE, as a dispersion-managed soliton (DMS). 
Importantly, while for the NLSE and the DMNLSE one can find solutions for any real A, for 
the CGLE and the DMGLE nontrivial solutions exist only for discrete values of X. 

Figure 2 shows the real and imaginary parts of stationary solutions of the NLSE, CGLE, 
DMNLSE and DMGLE, while figure 3 shows the modulus of the same solutions in semi- 
logarithmic scale. The values of A that yield the unique solution with the given value of the 
parameters for the CGLE and the DMGLE are respectively A^gie = 3.86 and A^mgie = 4.43, 
as obtained numerically using the methods described in appendix. For the NLSE and the 
DMNLSE, the value of A was chosen so that the corresponding solutions have the same peak 
amplitude as those of the CGLE, respectively. Such values are easily obtained by noting that, 
for both the NLSE and DMNLSE, A = maxxem. |m(x,0)|. In all cases, 1024 Fourier modes 
were used in the simulations. 

It should be clear from these figures that the solutions of the DMGLE combine some 
features of the solutions of the CGLE with some of the solutions of the DMNLSE. For 
example, it is evident from both figure 2 and figure 3 that the shape of the soliton solutions 
of the DMGLE is remarkably similar to that of the soliton solutions of the DMNLSE. At the 
same time, however, figure 2 shows that, like the solutions of the CGLE, and unlike those 
of the DMNLSE, a small but nonzero imaginary component appears. As a consequence, the 
solutions of the DGMLE do not appear to possess zeros, unlike those of the DMNLSE. 
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Figure 2. Shape of the solutions of (5.2). SoUd blue line: CGLE (s = 0, g = l,b = c = 0.2); 
dashed black line: NLSE (s = b = c = g = 0); thick solid red line: DMGLE (i = 4, g = 2, 
= c = 0.2); dot-dashed black hne: DMNLSE (i = 4, 5 = c = | = 0). In all cases d=\. The 
value of X for the NLSE and the DMNLSE was chosen so that their solutions have the same 
peak amplitude as those of the CGLE and DMGLE, respectively. The imaginary part of the 
solutions of both the NLSE and the DMNLSE is identically zero, while the real part is almost 
indistinguishable from that of solutions of the CGLE and DMGLE, respectively. 
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X 

Figure 3. Absolute value of the stationary solutions in semi-logarithmic scale for the same 
cases as in figure 2. (Line colors and styles are also the same.) Solutions of the DMGLE 
have exponentially decaying oscillating tails, like those of the DMNLSE, with the frequency 
of the oscillations increasing along the tails. Unlike the solutions of the DMNLSE, however, 
solutions of the DMGLE do not appear to posseses zeros, due to the presence of a non-zero 
iinaginary part. 
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Parameter dependence. To further explore the similarities and differences between solutions 
of the DMGLE, CGLE and DMNLSE, we next look at their parameter dependence. Even 
though the parameter X appears in identical way in the integral equation for the solutions of 
the DMNLSE and DMGLE, it nonetheless plays a very different role in the two cases. To see 
why, we briefly look at the case s = and compare solutions of the NLSE and the CGLE. 
When transformed back to the physical domain, the solutions of (5.2) are, in this case, 

Mnise(^,0 = s,Qch{ax)&'''^^'l^ , (5.3a) 
Mcgie(-x,f) = \/d[A sech(77x)] '+'^e'M , {53b) 
where a is an arbitrary real constant, while [11] 



V ca Va a 

with 



-3{l-bc) + J9il-bc)^ + 8{b + c)^ 
a = 2v-bil-v^), v = ^_ .(5.4) 

The functional form of (5.3a) and (5.3b) is similar. The key difference, however, is that 
(5.3a) represents a one-parameter family of solutions, since a [the soliton amplitude, which is 
inversely proportional to the pulse width and coincides with the soliton eigenvalue in (5.2) as 
well as with the scaling parameter in (4.2)] can take any real value. In contrast, all parameters 
in (5.3^) are uniquely determined by the coefficients in the CGLE. Thus, none of them can 
play the role of the scaUng parameter in (4.2). Instead, that role for (5.3^) is played by y/g. 

A similar situation arises for the DMNLSE and the DMGLE. Namely, for the DMNLSE 
the eigenvalue A is also the soliton amplitude and the scaling parameter, and is in one-to- 
one correspondence with the pulse energy (even though the relation is not simply of direct 
proportionality as for the NLSE [41]). For the DMGLE, in contrast, the value of A is 
completely determined by that of the coefficients in the DMGLE. Thus, in both the CGLE 
and the DMGLE, the addition of gain and loss dynamics seems to simply pick out one 
particular solution in the family of solutions of the corresponding conservative model (NLSE 
and DMNLSE, respectively), without significantly altering its form — at least for moderate 
values of the gain/loss parameters. 

The above observations are corroborated by looking at the relation between the energy 
of the pulse and its root-mean- square width, defined respectively as 



E = J\uf-6x, x=^fmjE, (5.5) 

where 

M2 = /x2|Mpdx. 

For the solutions (5.3a) and (5.3^) of the NLSE and the CGLE, these integral quantities take 
on the following values: 

Enise = 2Ja , •^ise = n/{V3a) , (5.6a) 

£cgle = 3JV|T//c, Tcgie = ;rv/T//3|, (5.6b) 

with rj uniquely determined from b and c via (5.4) as before. Eliminating a from (5.6a) 
and g from (5.6^) we then obtain the well-known relations between amphtude and width, 
respectively for the NLSE and the CGLE: 

£'nisefnise = 2nd/V3, ^cgiefcgie = V3ndr]/c. (5.7) 
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Figure 4. Root-mean-square (rms) width of the stationary solutions versus their energy for 
the NLSE. CGLE, DMNLSE and DMGLE (line colors and styles are as in figures 2 and 3). 
The data for the NLSE and the CGLE are obtained respectively from the first and the second 
of (5.7), while those for the DMNLSE and the DMGLE are obtained from numerical solution 
of the integral equation (5.2) for different values of X (DMNLSE) and g (DMGLE). The value 
of all the other parameters is the same as in figures 2 and 3. 



Equivalent relations were derived in [16] for the DMNLSE using a variational approximation 
with a Gaussian ansatz: 



where y — a/s. Eliminating a from (5.8) then yields the equivalent of (5.7). To the best of our 
knowledge, however, no variational approach has been developed for the CGLE that can be 
extended to the DMGLE. In this case, therefore, the relation between energy and width must 
be obtained by numerically solving the integral equation (5.2). 

Figure 4 shows the rms pulse width versus its energy for the NLSE, the CGLE, the 
DMNLSE and the DMGLE. We see that, as a result of the dispersion management, the effects 
of gain dynamics are reduced compared to the constant dispersion case, and the similarity 
between the parameter dependence of the solutions of the DMNLSE and the DMGLE is 
even closer than that among the solutions of the NLSE and CGLE. [Note also that, for the 
DMNLSE, the data from (5.8) is almost indistinguishable from those coming from numerical 
solutions of (5.2).] Note however that the curves for the NLSE/DMNLSE are obtained 
by varying A, while those for the CGLE/DMGLE by varying g. Indeed, it is crucial to 
realize that each point in the curves for the CGLE/DMGLE is a solution of a different partial 
differential equation. As a result the above comparison may appear to be somewhat artificial 
at first. In practice it is not, however. This is because solutions with different amplitude are 
experimentally produced by varying the pump power [57], which has precisely the effect of 
changing the value of the linear gain coefficient — and thus of ^ — in the DMGLE (3.7). 




2a {l-v/TT7ln[(l + v/rT7)/3']/2}3'' 
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Stability. In the absence of gain dynamics, solutions of the NLSE and one branch of 
solutions of the DMNLSE are stable under perturbations. In the presence of linear, band- 
limited gain and gain saturation, one would expect the pulse solutions of the CGLE/DMGLE 
to be more stable than those of the NLSE/DMNLSE. Indeed, solutions of the CGLE with 
constant coefficients are stable for a > [11]. (Such is the case for all the solutions in 
figures 2-4.) We expect that in these situations, amplitude perturbations will be damped in 
the DMGLE like they are in the CGLE, making the stationary pulses of the DMGLE stable. 
But of course this conjecture must be verified via careful analysis. One way to do so would be 
to perform a linearized stability analysis, that is, to look for the eigenvalues of the linearized 
DMGLE operator L[w, u] in (4.6) when u is the solitary- wave solution (5.1) obtained via (5.2), 
along the lines of what was done for the DMNLSE in Refs. [17, 38, 54, 61]. 

6. Conclusions 

In summary, we have derived the dispersion-managed Ginzburg-Landau equation (DMGLE) 
as the equation that governs the long-term dynamics of systems described by the CGLE 
with time-dependent coefficients. In particular we have shown how the DMGLE arises 
in Ti:sapphire femtosecond lasers, we discussed the properties of the equation and of 
its solutions. Since the DMGLE is (like the CGLE, NLSE and DMNLSE) a universal 
equation, however, we believe that it will prove to be a unified model for the description 
many other experimental realizations of femtosecond lasers in addition to Ti: sapphire, such 
as similariton lasers [30], lasers using all-normal dispersion fiber [20], and those using 
waveguide arrays [56]. 

From a mathematical point of view, perhaps the most important feature of the DMGLE 
is that it is amenable to analytical treatment. In fact, the results of this work open up many 
interesting theoretical questions, as well as some important practical ones. In particular, the 
following are fairly natural open questions: 

• Whether a proof of existence of stationary solutions is possible (as was done in Ref. [66] 
for the DMNLSE). 

• Obtaining the asymptotic behavior of the stationary solution as x — > ±oo (as was done in 
Refs. [42, 43] for the DMNLSE). 

• Whether there exist various branches of solutions as a function of d as in the DMNLSE 
(e.g., see Ref. [18] and references therein). 

• Whether there exist multi-pulse solutions (as in the CGLE [12] and the higher-order 
DMNLSE [4]). 

• The stability of all of these solutions (e.g., see Refs. [54, 61] and references therein for 
the DMNLSE and [1 1, 32] for the CGLE). 

• Whether there exist "true" dispersion-managed solitons corresponding to the stationary 
solutions of the DMGLE (3.7); that is, whether there exist time-periodic solutions in the 
original, unaveraged system (2.1). (Recall that for NLS-hDM, radiative losses beyond 
all orders [64].) 

• Even if true dispersion-managed solitons did not exist in the original, un-averaged 
system, it should still be possible to prove that solutions of (2.1) remain 0{ta) close 
to those of (3.7) up to times 0{l/ta) (as was done for the DMNLSE [67]). 

• Characterizing the large-s limit of the equation and of its solutions (as was done in 
Ref. [3] for the DMNLSE). 
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• Formulating "slow" ordinary differential equations for the long-term evolution of the 
pulse parameters which bypass the fast dynamics inside each map (as was done in 
Ref. [15, 54] for the DMNLSE). 

• Characterizing pulse colhsions and interactions. Because neither the CGLE nor the 
DMNLSE [2, 37] are integrable, one would expect that the DMGLE is not either, 
implying that pulse interactions will be inelastic. 

• More generally, since the CGLE with constant coefficients describes a remarkably rich 
variety of physical phenomena, including chaos [58], it will be interesting to use the 
DMGLE to see how these phenomena are affected by the presence of dispersion and 
nonhnearity management. 

Since the Hamiltonian formalism is lost, however, settling these issues might be significantly 
more complicated than for conservative systems such as the NLSE and DMNLSE. 

A more complicated equation of DMNLS type with gain and loss was also recently 
studied in [5] as a model for Ti:sapphire lasers. Also, a review of different mathematical 
approaches for the study of dispersion management in optical fibers can be found in Ref. [18]. 
All of the works cited in Ref. [18], however, deal with conservative systems. To the best of 
our knowledge, Ref. [5] and the present work are the first to generalize those tools in order to 
study dispersion-managed systems with significant gain dynamics. 

With regard to more practical issues pertaining to Ti: sapphire femtosecond lasers, we 
first note that a multi-dimensional version of the DMGLE such as the one presented in 
section 2 could be useful in order to take into account the transverse dynamics of the pulses 
in the cavity. In any case, the DMGLE can now be used to study the sensitivity of pulses 
with respect to perturbations, and especially quantum noise, using the linear modes and their 
adjoints to guide importance-sampled Monte Carlo simulations along the line of Refs. [41, 49- 
51]. (For an introduction to importance sampling tailored to the study of noise in hghtwave 
systems See Ref. [51].) Such a study will be instrumental to determine the true stability 
properties of pulses in these lasers and consequently to obtain the comb linewidth, which is 
an important step [7, 46] toward determining the ultimate accuracy of femtosecond lasers as 
optical atomic clocks. 
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Appendix: Numerical methods for the DMGLE 

Two relevant issues are: (i) efficient methods for numerical integration of the time-dependent 
DMGLE (3.7), and (ii) numerical methods to find dispersion-managed soUton solutions, i.e., 
solutions of the nonlinear integral equation (5.2). The first issue is rather straightforward: 
equation (3.7) can be integrated using essentially the same techniques as for the DMNLSE. 
These techniques include a method for the fast evaluation of the double integral. Since this 
method was discussed in detail in the appendix of Ref. [41], we do not repeat that discussion 
here, and we turn our attention instead to the nonlinear integral equation (5.2). 
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When b = c = g = Q (that is, for the DMNLSE), (5.2) can be efficiently solved 
numerically using Petviashvili's method (as discussed in the Appendix of Ref. [41]). Thus, 
here we only need to address the case when ^, c or ^ are nonzero. For simplicity, we first 
present the method for the CGLE (that is, when s = 0), even though in that case one can find 
stationary solutions analytically. This will allow us to explain the relevant ideas — which are 
the same for both the CGLE and the DMGLE — without some of the notational comphcations 
resulting from the presence of the kernel r{k'k") in (5.2). 

As with the DMNLSE, we solve the nonlinear integral equation (5.2) numerically by 
introducing an appropriate iteration scheme. Compared to the DMNLSE, however, two 
additional complications must be addressed: 

(i) An appropriate correction factor is necessary for convergence, since a standard Neumann 
iteration diverges. This issue is similar to that arising for the DMNLSE. In that case, the 
problem can be solved using Petviashvili's method. That method does not converge 
for the CGLE/DMGLE, however, due to the presence of both a real and an imaginary 
component in the solution. 

(ii) The propagation constant n is unknown a priori. That is, unUke the NLSE/DMNLSE, 
soliton solutions of the CGLE and DMGLE only exist for certain discrete values of 
the propagation constant, cf. (5.3^) and (5.4). Since these values are not known in 
advance, in addition to looking for the functions and Fj the solution method must 
simultaneously look for the values of for which (A.l) admits nontrivial solutions. 

We discuss both problems below. 

We seek a complex-valued function f{x) = /rc + ifim (with /re and /im both even 
functions of x) and a propagation constant n = such that F{k) := J\f] satisfies (5.2), 
which, when s = Q (that is, for the CGLE) reduces to 

(Ai + \ {d- ib)k^) f + igf = (1 + ic)f[ffr] . (A. 1) 

Decomposing (A. 1) into its real and imaginary parts yields the following two-component real 
system of nonUnear integral equations: 

(A-F/il)F-CN(F)=0, (A.2) 

where 



I is the 2x2 identity matrix and 

-»-{^$, ^^')^ -(J T)- 

(Of course /re and /im are also both even functions of k.) 

Let us first address the issue of the unknown propagation constant. Suppose that ilg is 
the exact propagation constant, and — jj.o+AiJ.. Equation (A.2) then yields 

A;uF = CN(F)-(A-|-Mol)F. (A.5) 

We can therefore can obtain A/z as 

(F,CN(F))-(F,AF) 
^~ WJ) ^ 
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Here the inner product of two real vector functions is 

(F,G) = /F(;t)-G(/t)d/fc. (A.7) 

Of course (A.6) contains the solution F, which is not known exactly during the iteration. 
Nonetheless, (A.6) provides a way to update our estimate of the correct eigenvalue. 

We now turn to the issue of convergence factors. This problem can be dealt with using 
the spectral renormalization method [8]. To implement this method, we start by noting that 
/I is not sign-definite, and even in those cases where is positive, some iterations might 
accidentally yield a negative estimate. In order to avoid problems with zero denominators, it 
is therefore convenient to add and subtract the term IF from equation (A.2), rewriting it as: 

(A+l)F = (l-;u)F + CN(F), (A.8) 

We then introduce the real convergence factor 7 and rescale the solution as F = yV. Noting 
that N(F) = 7^N(V), for the new field we obtain the system 

V = (A+l)-i[(l-M)V + r2CN(V)] , (A.9) 

which provides the basis for the iteration. An equations for the convergence factor is also 
needed, of course. This is obtained by using (A.8) with F = 7V and taking its inner product 
with V, which yields 

r2 = (V,(A + Ml)V)/(V,CN(V)). (A.10) 

Combining the two results, we can then define an iteration scheme as follows. At n = 0, 
choose an initial guess for fl^^ + ifj^ and and set yC) = 1 and V(°) = F(o). Then, at the 
(n + l)-st step of the iteration: 

1. Update the convergence factor 7^"^ and the current estimate of the eigenvalue ^^"^ 
using (A.IO) [with V = V(") and = 11^"^ and (A.6) [with F = 7(«)v("' and ji = n^"')]. 

2. Update the current estimate of the solution using (A.9) with the newly updated values 

At each iteration, the nonlinear terms can be computed using the fast Fourier transform (FFT), 
which reduces the computational cost of each step from from 0{N^) to 0{N^ logN), where N 
is the number of grid points or Fourier modes. 

The method to find soliton solutions of the DMGLE works exactly in the same way as 
that for the CGLE above. The only difference is that the term f\fff*] in the RHS of (A.l) is 
replaced by the double convolution integral P[f,f,f* , r], where 

piMi ,f;,r] = II A {k + k')f2{k + k")f; (k+k'+ ky(k'k"}dk'dk" . 

The nonlinear term N(F) is modified accordingly. Apart from these changes, the method 
is exactly the same as for the CGLE. Moreover, as in the DMNLSE, the resulting double 
convolution integrals can still be computed efficiently using FFTs. As this issue was explained 
in detail in the appendix of Ref. [41], we dot duplicate that discussion here. 
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